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We consider a model with anomaly-free Abelian gauge axial-vector symmetry, which is intended 
to mimic the standard electroweak gauge chiral SU(2)l X U(1)y theory. Within this model we 
demonstrate: (1) Strong Yukawa interactions between massless fermion fields and a massive scalar 
field carrying the axial charge generate dynamically the fermion and boson proper self-energies, 
which are ultraviolet-finite and chirally noninvariant. (2) Solutions of the underlying Schwinger- 
Dyson equations found numerically exhibit a huge amplification of the fermion mass ratios as a 
response to mild changes of the ratios of the Yukawa couplings. (3) The 'would-be' Nambu-Goldstone 
boson is a composite of both the fermion and scalar fields, and it gives rise to the mass of the axial- 
vector gauge boson. (4) Spontaneous breakdown of the gauge symmetry further manifests by mass 
splitting of the complex scalar and by new symmetry-breaking vertices, generated at one loop. In 
particular, we work out in detail the cubic vertex of the Abelian gauge boson. 

PACS numbers: 11.30.Qc 
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I. INTRODUCTION 

Use of the general principle of spontaneous breakdown 
of continuous symmetry for the SU(2)l x U(1)y invari- 
ant field theory description of electroweak phenomena is 
a necessity. Hard symmetry-breaking mass terms of both 
gauge and fermion fields ruin, either directly or by virtue 
of the loops, decent high-energy behavior of certain scat- 
tering amplitudes. 

In this respect the standard 'Higgs' realization of spon- 
taneous breakdown of the electroweak symmetry deserves 
admiration. At the expense of introducing a sector with 
a condensing scalar doublet it defines the operational 
(renormalizable, weakly coupled) description of virtually 
all electroweak phenomena so far explored. 

Because of theoretical drawbacks of the Higgs mecha- 
nism there are alternatives and/or generalizations thereof 
[l| (and references therein). In the near future the Large 
Hadron Collider at CERN will harshly test all of them. 
The aim of this paper is to add to the already existing 
list yet another realization of spontaneous breakdown of 
chiral and gauge symmetry. For clarity we illustrate it 
on physically and technically transparent Abelian pro- 
totype. Comparison with the standard Abelian Higgs 
mechanism can easily be made at any stage by heart. 

The basic idea Q is simple but subtle: As in the stan- 
dard Higgs realization we employ the complex scalar field 
but with an ordinary mass term of a complex scalar. 
Hence, in accordance with common wisdom there will be 
no spontaneous symmetry breakdown in the scalar sec- 
tor itself. Subtle is the nonperturbative self-consistence: 
The chiral symmetry breaking fermion proper self-energy 
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is both caused and causes the symmetry-breaking scalar 
proper self-energy, both solely due to a strong Yukawa 
interaction. 

The subsequent question of whether a massless gauge 
field introduced by gauging the U(l) axial symmetry de- 
scribes a massless or a massive spin-1 particle is dynam- 
ical and it was in general answered by Schwinger [3J: If, 
within a given dynamics, the gauge field polarization ten- 
sor develops a massless pole, its residue approximates the 
gauge boson mass squared [J] . Both in the Abelian Higgs 
model and in the model presented below the massless 
pole in the gauge field polarization tensor is due to the 
'would-be' Nambu-Goldstone (NG) boson of the underly- 
ing global symmetry. While in the Higgs mechanism the 
NG boson is pre-prepared in the elementary scalar field, 
in our case it is a composite object of both the fermion 
and the scalar fields. 

Finally, it is easy to realize that due to the dy- 
namical generation of the symmetry-breaking pieces in 
fermion and scalar field propagators there are new effec- 
tive symmetry-breaking vertices between the mass eigen- 
states of the fields. 

The paper is organised as follows: After introducing 
the model by its Lagrangian we discuss the possibility of 
the dynamical spontaneous breaking of the chiral sym- 
metry and introduce the appropriate formalism. Then 
we explore the consequences of the global symmetry: the 
existence of the 'would-be' NG boson and its interaction 
vertices with other particles. In this part we only repeat 
the main results of our preceding paper Q so we do not 
go into a detail. After this we present some consequences 
of the local symmetry: the mass of the gauge boson, in- 
duced by the bilinear coupling with the 'would-be' NG 
boson, and the existence of the effective trilinear coupling 
of the gauge bosons. In the end we present our numeri- 
cal results, especially the possibility of an arbitrarily high 
ratio of fermion masses with Yukawa coupling constants 
being of the same order of magnitude, together with the 
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description of our numerical procedure. 



II. THE MODEL 



In the full generality the model is defined by the La- 



grangian 
C 



1 



H~ Yukawa j 



w (D 
1 



(1) 



where the Yukawa interaction Lagrangian reads 

^Yukawa = V\ {^IL^XR^ + IpXR^XL^) 



(2) 



The Lagrangian possesses various symmetries. First, it is 
invariant under the global U(l)v 1 x U(l)v 2 x U(1)a sym- 
metry. The two global vector U(l) symmetries are gener- 
ated by the fermion number operators of the fermions 1 
and 2 and represent separate conservation of the number 
of both fermion species. The two naive independent axial 
U(l) transformations of the fermions 1 and 2 are related 
to each other by the Yukawa couplings to the scalar field 
<f>. Second, the Lagrangian is invariant under the local 
(gauge) U(1)a transformations, acting on the fields as 



-2i6(x) 



A„ + —d fJi 6(x) , 



(3) 



where the axial charges Q%, Q2 are defined as 

Qi = +1, 
Q2 = -1. 



(4) 



Considering two fermion species with opposite axial 
charges provides the simplest solution to an important 
theoretical constraint imposed on the model, which is 
the absence of the axial anomaly. The Yukawa coupling 
constants yx, yi are arbitrary real numbers. For the sake 
of simplicity we do not consider a more general Yukawa 
Lagrangian, allowing them to be complex. 



III. GLOBAL SYMMETRY 
A. Self-energies & masses 

In the first step we switch off the gauge interaction 
(g = 0) and demonstrate the spontaneous breakdown 
of the global U{\)a symmetry by finding symmetry- 
breaking proper self-energies in the fermion and scalar 
field propagators. We follow the general self-consistent 



nonperturbative method suggested by Nambu (and Jona- 
Lasinio) @,@|. 

First of all, it is convenient to introduce the Nambu- 
like doublet field for the scalar 



(5) 



$ = 

and its matrix propagator 

iD{x-y) = (0|T{*(x)$t(y)}|0) 



(6) 



The point is that this formalism allows us to treat the 
symmetry-preserving (the diagonal entries) as well as the 
symmetry-breaking (the off-diagonal entries) scalar prop- 
agators on the same footing. Similarly in the fermion 
sector, it is convenient to deal with the propagator 



iS(x-y) = (0\T{iP( x )$(y)}\0), 



(7) 



as it also contains both the symmetry-preserving (e.g. 
(4>l4>l)) and the symmetry-breaking (e.g. (iPl^r)) 
fermion propagators. 

However, although it would be possible, and even de- 
sirable, to deal with the propagators in the full general- 
ity, we will adopt for the present purposes the specific 
Ansatz: 



D-\p) 



p 2 -M 2 -LV) 
-U*(p 2 ) p 2 -M 2 



(8) 



Here the functions £j(p 2 ) and H(p 2 ) (in the follow- 
ing text often denoted simply as Sj >p and U p ) are the 
one particle irreducible (1PI) parts of the correspond- 
ing symmetry-breaking propagators, i.e. the proper self- 
energies. By this Ansatz we focus only on the symmetry- 
breaking parts of the propagators and totally neglect any 
possible effects of wave- function renormalization. For the 
purposes of the present paper this approximation is suf- 
ficient, since we wish only to demonstrate spontaneous 
symmetry breaking and not to make at this stage any 
phenomenological predictions. In a more realistic version 
of the model aspiring to describe the Nature and make 
reasonable physical predictions it would be, however, a 
must to include systematically all radiative corrections. 

Our general strategy in demonstrating the spontaneous 
breakdown of the U(1)a symmetry will be to search for 
the symmetry breaking parts of the propagators, i.e. the 
self-energies Ej and n. The key observation is, however, 
that at any finite order of perturbative expansion the 
U(1)a symmetry remains preserved and the self-energies 
vanish. The spontaneous symmetry breaking is a nonper- 
turbative effect. To treat it one has to employ some non- 
perturbative technique. We have chosen to make use of 
the Schwinger-Dyson equations, which represent a formal 
summation of all orders of perturbative expansion and as 
such they provide the desired nonperturbative treatment. 
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FIG. 1: The diagrammatical representation of Schwinger- 
Dyson equations ([9j|. The first line holds for both tpi and 
■02- The grey blobs stand for the proper self-energies, while 
the solid black blobs denote the full propagators. The double 
dashed line is for the Nambu $ doublet. 
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See also the diagrammatical representation of these 
Schwinger-Dyson equations - Fig. [T] If nonzero so- 
lutions £j and II exist, then they should necessarily 
be ultraviolet (UV)-finite since the corresponding counter 
terms are prohibited by symmetry. 

Fermion masses mj are determined in terms of Ej by 
solving 

m 2 =X 2 (p 2 = m 2 j ). (10) 

By dimensional arguments the solutions must have the 
form 

m j = Mf j (y u y 2 ), (11) 

to be compared with rrij — (—2M 2 /A) 1 / 2 yj in the case 
of a condensing <f>. This is the crucial point: while in the 
standard Higgs mechanism the functions fj are simple 
linear functions, in the present model we expect some 
more complicated functions, which should also be non- 
analytic, because we deal with nonperturbative effects, 
that can not be obtained at any finite order of perturba- 
tive expansion in coupling constants yj. Moreover there 
is a possibility that a small change in yj (say, about one 
order of magnitude) might produce much larger change 
in order of magnitude of fermion masses. 

Like in the fermionic sector, the scalar spectrum can 
be obtained as a pole of the propagator: 

Ml 2 =M 2 ±\U{p 2 =Ml 2 )\. (12) 



The system of Schwinger-Dyson equations forms an 
infinite tower of equations for all Green's functions of 
the theory. For practical calculation one usually has to 
truncate it at some level. We truncated it at the level of 
three-point Green's functions, which we approximate by 
the bare ones. The resulting system of Schwinger-Dyson 
equations in terms of the unknown functions Ej and II is 



(9) 
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This is easily interpreted: As a result of spontaneous 
symmetry breaking there appear two real scalar fields 
with different masses Mi. 2, instead of one complex field 
with mass M. 

Our numerical analysis confirms the existence of the 
nonzero solutions Ej and II. At present we are, how- 
ever, able to say nothing about their uniqueness. Upon 
performing the Wick rotation we have found numerically 
the real solutions Ej and IT with the following proper- 
ties: First, they vanish very fast at large p 2 (faster then 
any power of p 2 ). Second, the solutions are found so far 
only for large values of the Yukawa coupling constants yj 
(or, more precisely, for y\ and y 2 not being simultane- 
ously small). Third, fermion mass ratio m 2 /m 2 exhibits 
tremendous amplification upon decent changes of y\/y 2 . 
For example, for yi/y 2 — 77.4/88 we get m\/m 2 2 = 10~ 2 . 
Obviously this is alluring. If justified and understood 
analytically, this property of the solutions of the model 
should not be called fine tuning. Finding the explicit 
form of the functions fj (yi,y 2 ) in Eq. (jlip is our ulti- 
mate dream. 

More details on the numerics and adopted approxima- 
tions can be found in the Sec. [V] 



B. Nambu-Goldstone boson 

With the gauge interaction still switched off we reveal 
in the second step the Nambu-Goldstone (NG) excitation 
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FIG. 2: The diagrammatical representation of pole parts of 
proper vertices F^. and F^ (|17[) . The pole itself is interpreted 
as a propagator of intermediate massless scalar particle - the 
Nambu-Goldstone boson, depicted by solid double line. The 
small empty circles are vertices , P$ (|18(l . The external 
vector boson lines are indicated for the case the U(1)a sym- 
metry was gauged. 



and compute its effective couplings with fermion and bo- 
son fields. The analysis is standard: 

The conservation of the axial- vector current 

f A {x) = ^17^5^1 " ^27^75^2 + 2i [(d^)U - rfd^] 

(13) 

implies the axial-vector Ward identities for the proper 
vertex functions: 

(p + q,p) = Qj [Sj\p + q) 7 5 + J 5 Sj\p)] , 
q^(p + q,p) = -2D- 1 (p + q)~+2~D- 1 (p), 

(14) 

where the proper vertex functions T^, . , correspond to 
the Green's functions 



G» (x,y,z) = (0\T{j^(x)^(y)^(z)}\0) , 



(15) 



G${x,y,z) = (0|T{^(x)$(y)$t( z )}|0) 
with the full propagators of the external legs cut off. The 



1 

-1 



(16) 



operates in the <f> — 0' space and is quite analogous to 75 
in the fermion sector. 

For Tjj and II different from zero these identities im- 
ply that the proper vertices themselves develop poles for 
q 1 -0: 



r§, po ie(p + g,p) = ^[-2D- 1 (p + q)E + 2ED- 1 (p)] . 



(17) 



In fact the poles are nothing but the manifestation of the 
presence of the NG boson, see Fig. [2j These identities 
also allow to explicitly compute the interaction vertices 
of the NG boson with the fermions and scalars; for the 
details see Ref. [2[, the result is [within the Ansatz ©] 



N 



3 [£> + g) + S»] 75 , 



P*(p + q,p)= (18) 

_2_/ u(p + q) + n( P ) 

N \-U*(p + q)-Il*(p) 
The normalization factor N can be expressed as 



N= /^(0) + ^ 2 (0) + J (O), (19) 
where the functions J 1 p j and are defined by 



-iq^dq 2 ) 



d 4 fc (k - qY^ 3 ,k Sj.fc + ^j,k-g 
d 4 fc {2k-qY{k 2 - M 2 ) 



Re 



n* (n fc + n fe _ g ) 



(20) 



(2tt) 4 (fc 2 - M 2 ) 2 



Af 2 
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IV. LOCAL SYMMETRY 

A. Gauge boson mass 

In this step we switch on the gauge interaction pertur- 
batively and compute the axial-vector gauge boson mass 
squared as a residue at the massless pole of the gauge 
field polarization tensor Ii^ u (q). As in other strongly 
coupled models in which the NG boson is a composite 



[E S S B B we are only able to compute the longi- 
tudinal part of rP y due to the bilinear coupling of NG 
boson with gauge boson - see Fig. [3J 

Insisting on transversality of II M1/ due to the axial- 
vector current conservation, we conclude that the gauge 
boson mass squared is softly generated in terms of Sj 
and II as 

Mi= 5 2 [J^(0)+/ v , 2 (0)+^(0)] ) (21) 
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VWD CKAA^ 

FIG. 3: The effective bilinear coupling of NG boson with 
gauge boson, giving rise to the longitudinal polarisation of the 
latter and hence generating its mass. The bilinear coupling 
is induced by the couplings of NG boson with fermionic and 
scalar fields P^. and P^, compare with Fig. 
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(22) 



For q 2 = we have explicitly 



/*,(0) = -8i 
/*(0) = 16i 



d 4 fc S i , fc [E i , fc -2fc 2 a | T E,- fe ] 



(2tt)4 

d 4 fc fc 2 |n fc | [|n fc | - (fc 2 - 7v/ 2 )^|n fc |] 



fc2 - *h 



(2tt) 4 



(k 2 -M 2 ) 2 - |n fe | 2 



(23) 



In the weakly coupled Higgs model with the elemen- 
tary 'would-be' NG boson the spontaneous gauge boson 
mass generation is more transparent: The polarization 
tensor II'"' (g) can easily be computed completely and 
decent behavior of scattering amplitudes with longitudi- 
nally polarized massive gauge bosons can be illustrated 
explicitly. 

On the other hand, in strongly coupled models the lack 
of the bound-state spectrum apart from the NG sector 
guaranteed by the existence theorem of Goldstone im- 
plies: First, the g^ v part of the polarization tensor can 
hardly be computed. Second, decent behavior of scat- 
tering amplitudes with longitudinally polarized massive 
gauge bosons cannot be determined from first principles. 




FIG. 4: The full A vertex, denoted in the main text as 
iT MI/p (p, fc). Energy-momentum conservation is assumed: 
p + k + q = 0. 



The gauge interaction of A M with the massive fermions 
gives rise to the effective A 3 vertex depicted in Fig. [4] 
Up to g 3 the result (calculated as a proper vertex) can 
be written as the following sum (compare with Fig. 



B. A A vertex 

The spontaneous breakdown of the Abelian gauge chi- 
ral symmetry found in the low-momentum behavior of 
propagators should by construction manifest in nonin- 
variant loop-generated vertices. Being ultraviolet-finite 
they represent the genuine quantum-field theory predic- 
tions of the suggested nonperturbative approach. In the 
Abelian prototype presented here it suffices to provide 
one illustrative (still gedanken) example: 
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where (with the iO + 's suppressed) 



(24) 



-p-k)+iT^(k, 
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— Ct^ = Yl [ 
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+ 2 
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FIG. 5: Diagrammatical representation of the expression (1241) . 
Each diagram has the same assignment of momenta p, k, q 
and Lorentz indices n, v, p as the full vertex in Fig. [4] The 
symmetry factors are indicated explicitly. 



{Tr[7^7-(/-^) 7 ' ) (/+^ 75 ] + 



(2^) 4 [(^ 2 - M 2 ) 2 - |n,| 2 ] [ ((* + p) 2 - M 2 ) 2 - |IT, +P | 2 ] [ ((£ - k) 2 - M 2 f - \IU-k\ s 

[(f-M 2 ) 2 [n? +p n £ _ fe - c.c] + 
+ [(£ + pf - m 2 ] 2 [nj_ fc n< - c.c] + [(£ - k) 2 - m 2 ] 2 [n*n £+p - c.c] } , 



iTr p (p,k) = V 



d*£ 



g» u {2£ + p- k)P 



nj +p n^_*. - c.c. 



(2tt)4 [((i+ p y -M 2 f - |n /+p p] [((*- fc) 2 - m 2 ) 2 - 



(25) 



Being a combination of three axial vectors, the amplitude 
should change its sign under the operation of parity. In- 
deed it does: in the fermionic contributions it is due to 
the presence of the Levi-Civita tensor, while in the scalar 
contributions it is due to the following behavior of the 
scalar self-energy under parity: 



n Z ir 



(26) 



For illustration, let us now evaluate the amplitude 
iT^ vp (p,k) under certain approximations. First, we set 
the self-energies to be some constants. Note that the 
integrals remain perfectly UV finite. In fact, the scalar 
integrals vanish, while the quadratic divergences in the 
fermion integrals cancel each other since 



(27) 



(this, together with the requirement J2j Qj = 0; i s P re ~ 
cisely the condition for the axial anomaly to vanish) . The 



linear divergences cancel because of the symmetric inte- 
gration. Note that within this approximation the result 
is exactly the same as in the Higgs mechanism since we 
totally omit the nontrivial structure of the propagators, 
represented by nonconstant self-energies, and approxi- 
mate them by the pole contribution. Second, for the 
sake of simplicity, we put all external momenta on their 
mass-shell: 



p 2 = k 2 = q 2 = Mj 



(28) 



The energy-momentum conservation enables us to easily 
compute the dot products of external momenta: 

1 



p-q 



-.Mi 



(29) 



Constant fermion self-energies turn out to be the 
masses: Tij(p 2 ) — const. = mj. The resulting amplitude 
reads 



iT^ p {p, k) = G 



P e a/3 K 6 a/3 \P+ K ) 6 ap 



P 



(30) 
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or equivalently 



T^ p {p, k) = G \{q p k a - k»q a )p (3 e l ' p a0 + (p v q a - q" p 01 )^ e pp a0 + (k p p a - p"jfc Q )/e MI/ , 



a/3 



r 



(31) 



The latter form is perhaps more convenient because the 
invariance under exchange 



(32) 



is more apparent. The amplitude pip can be generated 
via the following effective Lagrangian: 

£ ff = Ge a ^ s {d a A a ){d?A°)(diA & ) . (33) 

Note that in the result l[3"Tj) there are no terms propor- 
tional to ^ p a: although they are present in the original 
\T^ vp [see ([25])]. This is a simple consequence of the re- 
striction p 2 = k 2 = q 2 , for in this case there is no way 



how to write down any nontrivial linear combination of 
terms e pi " p a p' 7 , t tiVp a k a and t tiVp a q a which would be in- 
variant under the exchange (|32p. For general momenta 
p 2 7^ k 2 ^ q 2 these terms would be present. 

The effective coupling constant G can be expressed as 



G = g 3 yQ 3 1 f(m 2 ,M 2 A ) 
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(34) 



The function / is defined by the integral 



/(m 2 ,M 2 



ir 2 M 2 



dz 



z(l-z) 



■. arctan 



j 



y/z(3z-4) + 



4m 2 



(35) 



r 



(here m 2 should be replaced by m 2 — i0 + whenever the 
correct branch choice of a multivalued analytic function 
is in question), which can be calculated analytically in 
some special cases: 

Ul " " 2A *] m2 ' (36) 
/(0,M 2 ) - 



6tt 2 M 2 ' 

More information about the shape of / can be extracted 
numerically - see Fig. [5] 



to be done. First, we switch to the Euclidean metric 
via the Wick rotation. By this we get rid of some poles 
in the fermionic and scalar propagators. The absence of 
poles makes numerical integration much easier. More- 
over, one can consider the self-energies to be real, with- 
out loss of generality. Note, however, that not all poles 
are removed, namely the pole in the scalar propagator 
remains: It takes place if H(p 2 ) = p 2 + M 2 for some p 2 . 



V. NUMERICS 

In this section the results of the numerical solution of 
Schwinger-Dyson equations ([9]) are presented. Moreover, 
as we consider the results quite interesting and impor- 
tant, the description of the numerical procedure together 
with discussion of its stability is presented. 



Second, we set A = 0. In fact, in the present model 
the A-term is not of high importance for the mass gen- 
eration. In the Higgs case the scalar self-interaction is 
of utmost importance: It fixes the numerical value of 
the condensate and stabilizes the scalar field energy. In 
our case both reasons for considering A ^ are clearly 
absent. However, in full treatment including the pertur- 
bative effects of renormalization yet to be done the term 
\\(<ft 4>) 2 will become indispensable as a counter term. 



Numerical results 



The system of equations © is quite difficult to be 
solved even numerically, so certain approximations have 



As a net result, we solve the following set of equations 
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FIG. 6: The independence of the function /(m 2 , M 2 ). Both 
quantities are normalized by m 2 to be dimensionless; note 
that m 2 /(m 2 ,M 2 ) is only function of M 2 /m 2 . The cusps 
appear at M 2 — 3m 2 and M 2 = 4m 2 . The former one also 
indicates the beginning of the imaginary part. 



FIG. 7: Typical shape of the solutions Ei(p 2 ), E2(p 2 ) and 
II(p 2 ) to the system of equations (|37jl . computed here for 
Yukawa coupling constants y\ — 83 and t)2 = 88. Note the 
saturation of the self-energies at low momenta and fast de- 
crease at high momenta. 



for the unknown functions £i(p 2 ), ^{p 2 ) and H(p 2 ): 
d 4 fc Ei fc Uk 



Vi 2 



d 4 fc £ 2 , fc 



n 



k—p 



(2tt) 4 fc 2 + E| fe [(fc - p )2 + M 2]2 _ H 2 fe _ p 

£j,fc-p 



d 4 & Sj,fe 



j=l,2 



(2^) 4 P + £ 2 fc (fc-p) 2 + £ 2 fc _ p 



(37) 



where p 2 = p 2 , + Pi + 2?2 + P§ > 0. 

The solutions can be classified with respect to the pairs 
of coupling constants (y%, 1)2) ■ The mass parameter M is 
not relevant, since as the only parameter with dimension 
of mass in the theory it serves just as a scale parameter 
for the self-energies and momenta. The typical shape of 
the resulting self-energies is depicted in Fig. \7\ They 
are saturated at low momenta and fall down rapidly at 
high momenta. Using logarithmic scale on y-axis in Fig. 
[7] it would be apparent that the self-energies fall down 
actually faster then any power of p 2 . 

We have found interesting qualitative differences be- 
tween the solutions of Eq. (|3"Tj) for different values 
(yi , 2/2)- The main results are summarised in Fig. [8J 
There is an area in the (t/i, 2/2) plane, located around the 
axis yi = 2/2 and denoted as (II), where both Si and £ 2 
are nonzero, and areas where only one of them, £1 or 
£2, is nonzero [areas (I) and (III)]. There is also an area 
(IV) where the pole in the scalar propagator, mentioned 
earlier, matters. The self-energies (which are expected to 
be imaginary here) were not computed in this area, we 
do not know anything about their behavior here. 

Our aim was to find the dependence of the spectrum - 
the masses of the fermions, scalar bosons and the vector 
boson - on the Yukawa coupling constants 2/1 and 2/2 [or, 
more precisely, on absolute values \yi\ and I2/2I, due to 



the shape of the Schwinger-Dyson equations [)37p]. For 
the calculation of masses we have used Eqs. (TTU1) . (|I2I) . 
(f2"Tj) . in the case of vector boson Wick rotated to the 
Euclidean metric. We have probed the j/i^-dependcncc 
along the cut depicted in Fig. [51 since it connects all the 
three main areas (I), (II) and (III) and therefore the re- 
sulting 2/1,2-dependence of the spectrum can be regarded 
as quite typical. The results are depicted in Fig. [9] and 
Fig. [TUJ Note how the critical lines between the areas 
are evident in the 2/1.2-dependence of the spectrum. The 
most significant result - the behavior of the fermionic 
spectrum - can be seen in Fig. 1101 As 2/1 approaches the 
critical line between (II) and (I) [or (III), respectively] 
in the direction from (II) to (I) [(III)], the ratio m\/m\ 
becomes arbitrarily high (low)! 



B. Numerical procedure 

Let us now say more about the numerical procedure 
itself. More formally, we have to solve the following set 
of equations: 



Ei = Fi[Ei,n], 
E 2 = F 2 [E 2 ,n], 
n = G[£i,E 2 ], 



(38) 



where Fj, G are some functionals. To solve this system, 
we have adopted the method of iterations. After choosing 
some initial Ansatz for E's: 



v (°) v(°) 
^1 i ^2 



(39) 



and calculating also the 'zeroth' iteration of the II: 

=G[E< 0) ,Ef], (40) 
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FIG. 8: The (3/1,3/2) plane with indicated areas of different 
behavior of the system of equations (|37p . According to the 
resulting fermion self-energies there are three main areas: first 
where Ei = and E2 7^ 0, second where Ei 7^ and E2 7^ 
and third where Ei 7^ and E2 = 0, denoted as (I), (II) and 
(III), respectively. There is also another important area, de- 
noted as (IV), where the pole in the scalar propagator takes 
place (and hence the imaginary parts of the self-energies are 
expected to appear). However we do not know anything about 
behavior of the self-energies in this area because of difficul- 
ties when numerically integrating the pole. The dashed line, 
going from 3/1 = 72 and 3/2 = 88 to y\ — 104 and 3/2 = 88, 
shows where the dependence of the spectrum on the Yukawa 
coupling constants was probed - see Fig. [5]and Fig. 1101 
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FIG. 9: The 3/1 -dependence of the scalar masses M lj2 with 
fixed 3/2 =88. 
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FIG. 10: The 3/1-dependence of the fermion and vector boson 
masses mf^ and M\ with fixed 3/2 =88. 



the iterating process is established: (n > 1) 
= Fxp^-^IT^- 1 )], 

s< n) = F 2 [z% l - 1 \n( n -V], 

n („) = gpMJ^M)]. 



(41) 



If this procedure converges, we have the solution. In 
order to control the convergence we used the quantity 



/(») = 



J* 



(n-1) 



(42) 



where X stands for any of £j, II. If I^ n > — > 1, we con- 
sidered XW to be the solution. This simple criterion is 
sufficient in our case because the functions X'"' turn out 
to be well behaved: there are no intersections between 
distinct iterations, or, loosely speaking, their shapes are 
similar, the only changes are in their 'size' (at least in 
the nonasymptotic region). 

Usual behavior of such a nonlinear system in the case of 
only one equation for one unknown function (for instance 
an equation for £ with II set to be a constant) is such that 



for any initial Ansatz the iteration procedure converges 
to some solution. The solution may be trivial (i.e. j' n J < 
1, I^ n > ■¥* 1) or nontrivial (i.e. I^ n > — > 1), depending 
whether y < y alt or y > y C rit, respectively, for some 
critical value y C rit- 

In our case of more coupled equations the situation is, 
however, different. Having some fixed Ansatz, there are 
three possibilities: (1) The iteration procedure converges 
to the trivial solution [which exists always for all (3/1, 2/2)]- 
(2) It blows up (i.e. JW > 1, -+> 1). (3) It converges 
to a nontrivial solution (if it exists), depending on the 
(2/1,3/2) (see Fig. El. 

This behavior clearly depends on the Ansatz. To be 
specific, our class of Ansatze consisted of functions 



if V) = 4° V) = 



(p 2 + M 2 ) 



2 \2 



(43) 



with x being a free parameter. For x too small or too 
large the iteration procedure goes to the trivial solution 
or blows up, respectively. Properly adjusting the value of 
x between these two extremes one can manage that the 
iteration procedure converges to the nontrivial solution. 
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C. Numerical stability 

A special care was taken to check whether the numeri- 
cal procedure is stable in the sense that its results (espe- 
cially the strong 7/ 12 -dependence of the fermion masses) 
remain unchanged upon the changes of the numerical al- 
gorithm. We tested four main variations of the algorithm: 

Class of Ansatze Besides (|43[) we considered also other 
Ansatz classes, all functions decreasing as some 
power of p 2 or exponentially, respectively. For some 
of them (some very rapidly decreasing exponen- 
tials) it was not possible to adjust the variable x 
[see (|4"3")) ] to converge to any but the trivial solu- 
tion. For the other Ansatz classes, however, the 
results were nontrivial and they all coincided. 

Step size There is, of course, a step size dependence 
of the numerical integration results, the important 
question is, however, how this dependence behaves 
for arbitrarily small step sizes. If there is no sen- 
sible (i.e. finite) limit of the integral as the step 
sizes are going to zero (the continuum limit), the 
results of the numerical integration have no mean- 
ing. We checked that this limit does exist and that 
all described phenomena, especially the strong yi,2- 
dependence of the fermion masses, are present in it. 

Integration method Having expressed the integral 
equations (|37|) in the hyperspherical coordinates, 
we can do two angular integrations analytically. 
The two remaining integrations (one over momen- 
tum and one angular) must be performed numeri- 
cally. For this purpose we employed consecutively 
the trapezoidal rule and the Simpson's rule, for 
the angular integration also the Gauss-Chebyshev 
quadrature formula (using Chebyshev polynomials 
of the second kind). The final results for all inte- 
gration methods agreed, the differences were only 
in the step size dependence, i.e. in the speed of 
convergence to the continuum limit. 

Momentum cut-off Since the momentum integral is 
over the infinite interval, for the purposes of the 
numerical integration a momentum cut-off must be 
introduced. However it turns out that one does not 
need to check the cut-off dependence of the results. 
That is because the resulting self-energies are so 
rapidly decreasing that the contribution from the 
high momenta is negligible. 

Moreover, in order to check the consistency of our nu- 
merical method by a comparison with an independent 
result, we calculated the equation for Ej [see (|37|)] with 
n set to be a constant [up to our knowledge, there are 
no independent calculations of the full set of the coupled 
equations (|3T[) we could compare with] and compared 
our result with the results of [12| [Eq. (2.11) and Fig. 3 
therein]. They coincided. 



VI. CONCLUDING REMARKS AND OUTLOOK 

Within the effective quantum field theory for elec- 
troweak interactions at forthcoming energy scale we find 
the use of scalar fields still welcome: Their Yukawa 
couplings with massless elementary fermion fields break 
explicitly the huge chiral symmetries of three elec- 
troweakly identical fermion families exactly to the sancti- 
fied SU(2)l x U(1)y- Alternatives without scalars strug- 
gle with guaranteeing unobservability of physical conse- 
quences of these unwanted symmetries. 

Linear dependences of the wildly wide, sparse and ir- 
regular fermion mass spectrum upon the Yukawa cou- 
pling constants of the Standard model are, however, a 
drawback. It is feasible that the fermion masses can be 
generated by the Yukawa interactions dynamically with- 
out ever referring to the scalar field condensation [7|, ll3| . 
In such a case the dependences of the fermion masses 
upon the Yukawa coupling constants would be nonana- 
lytic. In an Abelian prototype we have provided a bona 
fide indication that this is possible. Most importantly, we 
have convincingly demonstrated that a small change in 
the ratio of the Yukawa couplings may lead to an orders- 
of-magnitude change in the ratio of the generated fermion 
masses. We hope that this mechanism could provide a 
natural explanation of the hierarchy of fermion masses in 
the Standard model. 

The price we pay is that the Yukawa interactions hav- 
ing the expected properties must be strong. But we 
are accustomed to live with QCD strongly interacting 
at large distances even though we still don't know how 
to solve it (incidentally, apart from the NG sector). We 
are tempted to interpret our findings as manifestations 
of nontrivial fixed points [l4j in our model. 

The bonus we get is an interesting relation of fermion 
and gauge boson mass generation mechanisms. Indeed, 
in contrast with the Higgs mechanism in our case for the 
Yukawa couplings set to zero both fermion and gauge 
boson masses vanish. 

We are confident that with some technical modifica- 
tions the generalization of the Abelian mechanism pre- 
sented above to SU(2)l X U(l)y can be done [la ]. 
Whether it can be converted into a phenomcnologically 
viable alternative of the Standard model remains, how- 
ever, to be seen. 
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